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ABSTRACT 

We report on a 70ks XMM-Newton ToO observation of the newly discovered ac- 
creting millisecond pulsar, IGR J17511-3057. Pulsations at 244.8339512(1) Hz are 
observed throughout the outburst with an RMS pulsed fraction of 14.4(3)%. Pulsa- 
tions have been used to derive a precise solution for the Porb = 12487. 51(2)s binary 
system. The measured mass function indicates a main sequence companion star with 
a mass between 0.15 and 0.44 Mq. 

The XMM-Newton 0.5-11 keV spectrum of IGR J17511-3057 can be modelled by 
at least three components, which we interpret, from the softest to the hardest, as mul- 
ticoloured disc emission, thermal emission from the neutron star surface and thermal 
Comptonization emission. Spectral fit of the XMM-Newton data and of the RXTE 
data, taken in a simultaneous temporal window, well constrain the Comptonization 
parameters: the electron temperature, kTe=5lt4 keV, is rather high, while the optical 

depth (T=1.34to:o6) is moderate. 

The energy dependence of the pulsed fraction supports the interpretation of the 
cooler thermal component as coming from the accretion disc, and indicates that the 
Comptonizing plasma surrounds the hot spots on the neutron star surface, which in 
turn provide the seed photons. Signatures of reflection, such as a broadened iron Ka 
emission line and a Compton hump at ~ 30 keV, are also detected. We derive from 
the smearing of the reflection component an inner disc radius of ^ 40 km for a 1.4 
M0 neutron star, and an inclination between 38° and 68°. 

XMM-Newton also observed two type-I X-ray bursts, whose fluence and recurrence 
time suggest that the bursts are ignited in a nearly pure helium environment. No 
photospheric radius expansion is observed, thus leading to an upper limit on the 
distance to the source of 10 kpc. A lower limit of 6.5 kpc can be also set if it is assumed 
that emission during the decaying part of the burst involves the whole neutron star 
surface. Pulsations are observed during the burst decay with an amplitude similar to 
the persistent emission. They are also compatible with being phase locked to pre-burst 
pulsations, suggesting that the location on the neutron star surface where they are 
formed does not change much during bursts. 

Key words: stars: pulsars: individual: IGR J17511-3057, X-rays: binaries 



1 INTRODUCTION 
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The discovery of 244.8 Hz pulsations in the X-ray emis- 
sion of the newly discovered tr ansient source, IGR J17511- 
3057 l|Markwardt et all 120091 '). has brought to thirteen 
the number of Accreting Millisecond Pulsars (hereafter 
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Figure 1. 2-60 keV light curve of the IGR J17511-3057 outburst 
as observed by the PCU2 of the PGA aboard RXTE. Circles and 
crosses refer to observations pointed towards IGR J 175 11-3057 
(Obsid P94041) and XTE J1751-305 (Obsid P94042), respec- 
tively. As the two PCA fields of view greatly overlap, the de- 
tected pulse frequency is used to discriminate from which source 
the emission observed by the PCA was actually coming from (see 
text for details). To which source the emission is attributed is 
indicated at the top of the figure. The time at which the XMM- 
Newton observation was performed is also shown. 



of similar broadened lines from AMSPs has a peculiar im- 
portance in the context of NS accretors. It may constrain in 
fact where the magnetosphere breaks off the accretion disc 
flow and lifts off matter to the accretion columns. 

The spectral and timing capabilities of XMM-Newton 
already proved fundamental to investigate the X-ray emis- 
sion and pulse properties of AMSPs especially at low en- 
ergies, as well as to detect iron emission lines at a high 
statistics. To these ends, we have obtained a 70ks Target 
of Opportunity (ToO) observation of IGR J17511-3057, per- 
formed roughly one week after the outburst onset. We report 
on the spin and orbital properties of the source in Sec. 13.11 
and on its spectrum in Sec l3.3l We also show how it is possi- 
ble to take advantage of a simultaneous Rossi X-ray Timing 
Explorer (RXTE) observation, in order to constrain its high 
energy emission (Sec. 13.3]) . 

Also two type I X-ray bursts were caught by XMM- 
Newton, making this observation particularly rich (Secl4]). 
We observed burst oscillations during the bursts, at the 
same frequency of the non-burst emission. It is of key im- 
portance to observe burst oscillations from sources, the ro- 
tational state of which is known at a great accuracy. The 
comparison with the properties of the non-burst pulsations 
may in fact assess the physical mechanism that produces os- 
cillations during bursts. These have already been observed 
to be strongly phase locked to the non-burst oscillations, 
possibly suggesting that the simil arities in their for mation 
are more profound than expected (jWatts et al.ll2008l '). 



AMSP) discovered so far. Such high spin frequencies are 
believed to ensue from a long phase ( Si 10* yr) of ac- 
cretion of mass and angular momentum onto the neutron 
star (NS), according to the so-called recycling scenario 
IJBhattacharva fc van den Heuvellll99ll ). 

In this paper we mainly focus on the spectral proper- 
ties of this source. The broadband continuum of AMSP has 
been invariably found to be dominated by a power-law like 
hard emission with a cut off at an energy between 30 and 60 
keV . This emission is interp reted as Compto nization of soft 
photons in a hot plasma (see |Poutanenll2006l . and references 
therein, for a review of the spectral properties of AMSP). 
At energies of '^l keV two soft compone nts are also gener- 
ally found IIG icrlihski & Poutancn 2005, GP05 hereinafter, 
IPapitto et"aill2009 : Patruno et al. 2009,:) . The cooler one is 
attributed to the accretion disc emission, while the hotter is 
interpreted as thermal emission of the hotspots on the NS 
surface. In this scenario, the hotspot thermal emission con- 
stitutes the bulk of soft photons that Compton down-scatter 
the hot plasma. The X-ray spectrum of IGR J17511-3057 is 
found consistent with this phenomenological description and 
the outlined physical interpretation fits well the observed 
spectral properties. 

A couple of AMSP have also showed the typical 
clues of disc reflection, a Compton hump at '--^ 30 
keV (IGierliriski et ahl |2002|) and a Kg em i ssion iron line 
JPapitto et al.l l2009l : ICackett et all l2009l : IPatruno et~aD 
I2OO9I ). In particular the line observed from SAX J1808.4- 
3658 is compatible with the typically broadened shape 
caused by the relativistic Keplerian motion of the inner rings 
of the accretion disc, i mmersed in the deep gravitational well 
of the compact object (jFabian et al.ll 19891 ). The observation 



2 OBSERVATIONS AND DATA ANALYSIS 

The X-ray transient IGR J17511-3057 has been discovered 
by INTEGRAL on 200 9 September 21 140, d uring its Galac- 
tic Bulge monitoring (JBaldovin et al.l l2009l ) . At that time 
also the RXTE Proportional Counter Array (PCA) detected 
a rising X-ray flux activity from the Galactic Bulge. This 
emission was tentatively attributed to two known nearby 
transien t sources, XTE J1751- 305, an AMSP spinning at 
435 Hz ()Markwardt et al.ll2002l ). and GRS 1747-312. RXTE 
therefore pointed at the position of XTE J1751-305 to 
discover that the emission instead originated by a source 
spinning at ~245 Hz. T his made IGR J17511 -3057 the 
twelfth AMSP discovered (JMarkwardt et al-llioogl ). A foUow 
up Swift ToO observation constrained the source position at 
a few arcsec accuracy, obtaining a value ~20 arcmin away 
from XTE J1751-305, thus c onfirming the discovery of a 
new source IjBozzo et al.ll20ld ). A Chandra observat ion fur- 
ther refined the source position (JNowak et aLll2009l ). giving 
the most accurate available estimate as RA=17h 51' 08.66", 
DEC=- 30° 57' 41.0" (1 sigma error of 0.6 arcsec), which is 
the one we consider in the following. 



2.1 XMM-Newton 

XMM-Newton observed IGR J17511-3057 as a ToO obser- 
vation for 70ks starting on 2009 September 21.037, ~ 9 days 
after the source discovery. The XMM-Newton observation 
is indicated by a horizontal bar in FigllJ where the RXTE 
PCA light curve of the outburst is plotted. 

The EPIC-pn camera operated in timing mode to al- 
low the temporal resolution (30 ^s) needed to study the 



X-ray spectrum of the newly discovered AMSP IGR J17511-3057 3 



CD 
T3 

CD 

E 

re 
> 

CO 

CD 
_cn 

Q. 



300 



200 



100 




-100 



-200 - \ / - -2 



E 

U) 
CO 
3 
"D 

en 
CD 

a: 



-300 

50 100 150 200 250 300 350 

Mean orbital longitude (degrees) 

Figure 2. Pulse phase time delays induced by the orbital motion 
as a function of the mean orbital longitude. Longitude 0° refers 
to the passage of the pulsar at the ascending node of the orbit. 
Residuals with respect to the best fitting orbit model are also 
over-plotted at a magnified scale (rightmost y-axis) . The reduced 
chi squared of the model is 1.15 over 684 degrees of freedom. 



millisecond variability of the source, and to avoid spectral 
deformations due to pile-up. Also EPIC-MOS2 camera ob- 
served in timing mode, while EPIC-MOSl operated in Small 
Window to provide an image of the source. No evident back- 
ground due to soft proton flares was recorded during the 
observation, and the full exposure window is retained for 
scientiflc purposes. XMM-Newton detected two type I X-ray 
bursts during its pointing. When analysing the persistent 
emissiorij, we discarded 20s prior and 110s after the burst 
onset. This choice is fairly conservative as the e-folding fac- 
tor of the burst decay is r ~ lis (see Sec|3]). 

XMM-Newton data have been extracted and reduced 
using SAS v. 9.0. A concatenated and calibrated event file 
has been created with the task epproc, specifying the 
Chandra coordinates of the source. The EPIC-pn spectrum 
has been extracted from a 13 pixels wide region around 
the source position (RAWX=36), equivalent to 53.3 arcsec 
(which should encircle more than 90 per cent of the source 
counts up to 9 keVo) A similar stripe centred on RAWX=10 
has been used to estimate the background. The background 
subtracted 1.4-11 ke\^ persistent count rate recorded by the 
EPIC-pn is 43.9 c/s on the average. Pile up is therefore not 
an issue for EPIC-pn data. Spectral channels have been re- 
binned to have at least three channels per resolution element 
and 25 counts per channel. 

EPIC-MOSl data have been extracted using the task 
emproc to provide a 101" x 101" image of the source. Using 



^ Even if the source is an X-ray transient, we refer to persistent 
emission as the main body of the outburst emission, and to burst 
emission when the source exhibits thermonuclear flashes, so called 
type I X-ray bursts. 

^ See XMM-Newton Users handbook, issue 2.6, available at 
http:// xmm.esac.esa.int /external /xmm_user .support. 

see Sec l3.3l for details about the energy interval chosen for the 
purposes of a spectral analysis. 



the task edetect_chain, we estimate the source position as 
RA=17h 51m 08.55s, DEC=-30° 57' 41.7", with an uncer- 
tainty of 4.3 arcsec. This is consistent with the more ac- 
curate Chandra determination. EPIC-MOS 2 data have a 
significantly lower statistics than that of EPIC-pn and are 
therefore discarded from the spectral analysis. 

Refiection Grating Spectrometers (RGS) operated in 
standard spectroscopy mode and data have been reduced 
using the task rgsproc. We consider only first order spectra, 
re-binned in order to have at least 25 counts per channel. 
The average persistent count rate is 0.4 and 0.5 c/s for RGSl 
and RGS2, respectively. 

All spectral fits were performed using XSPEC v. 12.5.1 

2.2 RXTE 

Similarly to the other AMSPs, only one third of the X- 
ray flux of IGR J17511-3057 is emitted in the energy range 
covered by XMM-Newton. In order to better constrain the 
broadband X-ray emission, we take advantage of two RXTE 
observations performed simultaneously to the XMM-Newton 
pointing, and amounting to an exposure of 8.8ks (Obs. 
94041-01-02-08 and Obs. 94041-01-02-18). 

RXTE observations have also been used to ex- 
tract a light curve of the whole outburst. Starting on 
2009 September 13.849, RXTE observed IGR J17511- 
3057 for ~ 24.5 d, achieving a total exposure of 611 ks 
(~30% of the whole outburst; Obsid P94041). The 2-60 
keV light cu rve extrac ted from the PCU2 of the PGA 
(IBradt et al.M l993: Jah oda et al.l :2006'1 data is plotted in 
FigUl Background has been subtracted using the faint model 
(pca_bkgd_cmfaintl7_eMv20051128.mdl), regardless during 
the outburst the source crosses the threshold of 40 c/s/PCU 
above which a bright model should be used. Such a choice 
was made to avoid the unphysical flux discontinuity that 
would take place when switching from one model to the 
other. In Fig[T] filled circles represent observations pointed 
to the IGR J17511-3057 position (ObsId P94041), while 
crosses refer to pointings in the direction of the nearby 
source, XTE J1751-305 (ObsId P94042). These two sources 
are 19.67 arcmin away, and there is an overlap of nearly 
80% between the respective PCA pointings (the FWHM of 
the PCA collimator is '^ 1°). RXTE observations are sub- 
ject to the contribution of both sources, no matter which 
the instrument is actually pointed to. The detected spin fre- 
quency is therefore used to discriminate which source emits 
the observed fiux. As during the first observation reported in 
Figfflthe 245 Hz periodicity is detected, it can be attributed 
to IGR J17511-3057 even if the instrument was pointing 
towards XTE J1751-305. The overlap between RXTE ob- 
servations of these two sources becomes even more evi- 
dent when looking at the late part of the outburst. Around 
2009 Oct 8 (Day 26 according to the scale used in Fig[l] 
where time is reported in days since 2009 Sep 12.0), the 
X-ray fiux shows an increase of more than 70% with re- 
spect to the previous emission. Most importantly, 435 Hz 
pulsations are detected, clearly indicating a renewed activ- 
ity from XTE J 175 1-305, as it has also b een confirmed by 
a nar row pointed Swift XRT observation (jMarkwardt et al.l 
[2OO9I). Observations performed between 2009 Oct 8 and Oct 
10 can be therefore safely attributed to a dim outburst of 
XTE J1751-305, similar to those already displayed in 2005 
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Table 1. Best fitting orbital and spin parameters of IGR J17511- 
3057 as observed by XMM-Newton. Errors in parentiieses are eval- 
uated at Icr confidence level, while upper limits are given at a 3cr 
level. 



a sin j/c (It-s) 

Porb (s) 

T* (MJD) 

e 

f(Mi,M2,i) (Mq) 



0.275196(4) 

12487.51(2) 

55094.9695351(7) 

< 8 X 10^5 

1.07086(5) X 10-3 



y (Hz) 
i> (Hz/s) 



244.8339512(1) 
< 3 X 10-" 



Marc h IJGrebenev et al.ll2005l ') and 2007 April (JFalanga et all 
[2003). 

The fact that IGR J17511-3057 is in the Galactic Bulge 
not only forces to deal with a crowded field, but also with 
the Galactic Ridge emission. It is in fact easy to see from 
FiglU that the PCU2 count-rate stays at a level of ~ 6.5 
c/s even after the XTE J1751-305 activity episode is over. 
Emission at the same flux level was observed also at the 
end of the 2002 outburst of XTE J1751-305 (GP05) and 
owes to the Galactic ridge. This emission introduces an ad- 
ditional background to all the PGA observations. This is 
clearly indicated by the prominent 6.6 keV iron line that 
appears in PCA spec tra when this addition al background 
is not accounted for (Markwardt ct al.l l2009l ). As we show 
in Sec lS.SI to subtract this background is of key importance 
to ensure agreement between the EPIC-pn and RXTE-PGA 
spectral data. To estimate a reliable spectrum of this addi- 
tional background we therefore summed up all RXTE data 
from 2009 Oct 11.400 to 22.745 (for an exposure of 87.8 ks). 
In the following we refer to this additional background as 
the Galactic ridge emission. 

In order to extract a PGA spectrum, we have consid- 
ered only data taken by the PCU2, that is the only always 
switched on during the .RJf'T'iJ observations overlapping with 
the XMM-Newton pointing. Data processing has also been 
restricted to its top xenon layer, as it is the less affected 
by the instrumental background. Only photons in the 3-50 
keV band have been retained as this band is best calibrated. 
Data take n by the High Energy X-ray Timing Experiment 
(HEXTE. [Rothschild etHI 1 19981 ) in the 35-200 keV band 
has been extracted considering the cluster B, which was the 
only one to perform rocking to estimate the background, at 
the time the observation took place. 



3 THE "PERSISTENT" EMISSION 

3.1 The pulse profile 

In this paper the pulsations shown by IGR J17511-3057 
are analysed using only the XMM-Newton data. A tempo- 
ral analysis based on RXTE data is instead included in a 
companion paper (Riggio et al., 2010, submitted). In or- 
der to perform a timing analysis we retain all the 0.3-12 
keV energy interval covered by the EPIG-pn, and report 
the photons arrival time to the Solar System barycentre us- 
ing the S AS task haryc e n, con sidering the Chandra position 
given bv iNowak et al.l l|2009l ). The 244.8 Hz periodicity is 
clearly detected throughout the XMM-Newton observation. 



o 
O 




Figure 3. Pulse profile obtained by folding the 0.3-12keV EPIC- 
pn time scries in 64 bins around the best spin frequency estimate, 
u = 244.8339512, after having corrected data for the orbital mo- 
tion. The solid line is the best fitting five harmonics function 
described in text. The first (dashed line), second (dotted line) 
and third (dash-dotted) harmonic are also over-plotted. An offset 
equal to the average count rate have been added to the harmon- 
ics count-rate for graphical purposes. Two cycles are shown for 
clarity. 



To estimate the spin and orbital parameters we focus on the 
persistent emission. We first fold 100 s long int ervals around 
the be st estimate of the spin frequency given bv lRiggio et al.l 
(J2009f ). The pulse profiles thus obtained are modelled with 
a sum of harmonic functions: 



F{<f>) = Ao{l + 



k 



AkCos\k( 



(1) 



where ^o is the average countrate, = 2-nh'{t — To) is the ro- 
tational phase, and Ak and (j)^ are the fractional amplitude 
and phase of the k-th harmonic, respectively. Two harmon- 
ics are enough to model pulse profiles obtained over 100s 
intervals. The relevant orbital parameters (the semi-major 
axis of the NS orbit, x = asini/c, the orbital period, Porb, 
the epoch of passage at the ascending node of the orbit, 
T* , and the eccentricity e) are estimated fitting the time 
delays affecting the first harmonic phas es with respect to a 
constant frequency model (see Eq.(3) of lPapitto et al.ll2007l . 
and references therein). The observed time delays, together 
with the residuals with respect to the best fitting model, are 
plotted in Fig[2] while the orbital parameters we obtain are 
listed in Table [T] They are perfectly compatible with those 
estimated by Riggio et al. (2010, submitted) considering the 
RXTE coverage of the outburst, even if the uncertainties 
that affect them are naturally larger because of the shorter 
exposure they are calculated over. 

After the photon arrival times have been corrected for 
the orbital motion of the pulsar, we folded data in 16 phase 
bins over longer time intervals (1000s), and derived more 
stringent measures of the phases and of their temporal evo- 
lution. The statistical errors on the phases obtained fitting 
the profiles so obtained are summed in quadrature with 
the uncertainty introduced by the errors affecting the or- 
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Figure 4. Variation with energy of the fractional amplitude of 
the first (top), second (middle) and third (bottom) harmonic. 



bital parameters (see Eq.(4) in lPapitto et al.ll2007l ). As there 
are no significant differences between the estimates we ob- 
tain with the two harmonic components, we present here 
only the resuhs based on the fundamental one. The phase 
evolution can be successfully modelled by a constant fre- 
quency iXred — 1-02 for 68 d.o.f.), whose estimate is given 
in Table [T] The quoted error is evaluated considering also 
the uncertainty int roduced by the positional error-box (see 
iBurderi et al.l 120071 ). The introduction of a quadratic term, 
possibly reflecting a spin evolution of the source, is not sig- 
nificant, and it is possible to derive a 3a upper limit on the 
spin frequency derivative of 1 x 10~^^ Hz/s. 

Folding over the entire observation length the 0.3-12 
keV EPIC-pn time series we obtain the average pulse profile 
plotted in Fig|3l Five harmonics are needed to successfully 
fit (i.e. x^/dof=53.1/54) the pulse profile using Eq.(IT|), while 
AMSPs pulse profiles are generally fitted by just two com- 
ponentcl. Such a complexity is unveiled probably thanks 
to the high statistics and low background granted by the 
EPIC-pn, and, most importantly, to the larger pulse fraction 
IGR J17511-3057 shows with respect to the other objects 
of this class. Accounting for the average background count- 
rate (2.225(6) in the 0.3-12 keV band), the fractional ampli- 
tudes of these harmonic components are in fact, 20.21(8)%, 
2.66(8)%, 0.98(8)%, 0.39(8)% and 0.33(8)%, respectively, 
where the errors in parentheses are quoted at la confidence 
level. The total RMS pulsed fraction can be estimated as 
(Efc^'/2)'/' = 14.4(3)% RMS. 

The pulse of IGR J17511-3057 also shows strong spec- 
tral variability. The amplitude of the first harmonic increases 
with energy until it reaches an approximately constant value 
of ~ 0.22 at ~ 3keV. The second and third harmonic are in- 
stead more regular (see Fig. |4]). We show in Sec l5.3l how the 
decrease of the pulsed fraction can be understood in terms 
of the shape of the various components used to model the 



* We note that more than two harmonics has bee n detected dur- 
ing a subset of observations of SAX J1808. 4-3658 IjHartman et al.l 
I2OO8I) and XTE J1807-294 l lPatruno et al.ll2009l) 




Energy (keV) 

Figure 5. Phases in /is of the first (solid line) and second har- 
monic (dotted line) with respect to energy. An common offset has 
been set on phases so that the phase of the fundamental between 
1 and 2 keV is equal to 0. 



X-ray spectrum. Phase lags are also observed (see Fig[S}. In 
particular the phase of the fundamental shows an excursion 
of '~^ 230^s (0.06 in phase units) between 1 and 10 keV, 
with the pulses at low energy lagging those at higher ener- 
gies. Lags are also shown by the second harmonic, even if 
their significance is lowered by the larger error bars affect- 
ing these estimates, with respect to those calculated on the 
first harmonic. The behaviour with energy of the second har- 
monic phases is anyway more regular. While at low energies 
the second harmonic phase anticipates the first harmonic, at 
higher energies they become comparable. We thus conclude 
that the overall pulse shape is energy dependent. 



3.2 The XMM-Newton spectrum 

We first analyse the XMM-Newton spectrum, considering 
data taken by the two RGS (0.5-2.0 keV) and by the EPIC- 
pn (1.4-11.0 keV). EPIC-pn data below 1.4 keV are dis- 
carded as they show a clear soft excess with respect to RGS 
data, regardless of the model used to fi t the spectrum. Suc h 



an excess was alr e ady noticed by, e. g., Boirin et al.l 



laria et al 



D'Ai et al. 



20051): 



20091 ): 



2009 [ ): lD^i et al.l l|2009l ): lPapitto et al.l 1 
(2010;), analysing observations performed by the 
EPIC-pn in timing mode. The absence of such a feature in 
RGS data indicates how it is probably of calibration origin, 
hence the choice to retain only data taken by the EPIC-on in 
the 1.4-11.0 keV energy interval, for the purposes of a spec- 
tral analysis. The relative normalisation of the RGSl and 
RGS2 with respect to the EPIC-pn camera are left free. We 
find RGS1/PN=1.02(1) and RGS2/PN=0.98(1), regardless 
of the particular spectral model considered. The numbers 
quoted in parentheses are the errors at 90% confidence level, 
level at which the uncertainty on any spectral parameter is 
quoted in this paper. 

The 0.5-11 keV spectrum of IGR J17511-3057 is ener- 
getically dominated by a power law, which is easy to inter- 
pret as the realization of a broader Comptonized emission 
on a limited bandwidth. The first Comptonization model we 
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Table 2. Fitting parameters of the 0.5-11.0 keV combined RGS + EPIC-pn spectrum of IGR J17511-3057. 
Models are defined as follows: Model A, phabsX6dge(diskbb+bbodyrad+nthcomp) with \i.Taoft =kTBfl; Model B, 
phabsxedge(diskbb+bbodyrad+nthcomp) with kTj,o/t T^kT^a; Model C, same as model Model B with the addition of a diskllne 
to model the Ka iron emission feature. The disc outer radius is fixed to 1 X 10^ Rg. All these models are evaluated for an electron 
temperature of the Comptonized component fixed to 100 keV. Errors on each parameter are quoted at the 90% confidence level, as 
for all the spectral parameters in the rest of the paper. 



Model 


A 


B 


C 


nH (10^2 cm-2) 


1 in+o.07 


0.98 ± 0.04 


O.96IH2 


T O VIII 


0.22 ±0.08 


0.21 ±0.07 


0.20 ±0.07 


kTi„ (keV) 


p,r,i+0.02 


0.33J:H^ 


r. 34+0.04 


Ri„(cosJ)l/2 (dskm) 


63+28 


20+^ 


17+5 


kTBB (keV) 


n 30+0-02 


fi2+°°3 
U.0^_0 Q2 


m+°-°'^ 


Rss (dg km) 


7±1 


71+0.5 
'■^-0.6 


6.7+g-^ 


r 


1.710 ±0.008 


2.4 ±0.2 


2-llg-? 


kT,,jt (keV) 


=kTss 


1.5 ±0.01 


1-36+r? 


EFe (keV) 


— 


— 


"■'^-0.09 


/3 


— 


— 


4 1+0.8 


R„ {GM/c'^) 


— 


— 


371^1 


i(°) 


— 


— 


>37 


EW (eV) 


— 


— 


43.9 ±0.06 


2 


2008.1/1598 


1775.8/1597 


1740.7/1592 



consider is nthcomp IjZdziarski et al.lll996l : IZvcki et al.lll999l ). 
It describes a thermal Comptonization spectrum in terms of 
the temperature of the input soft photons (kTso/t), the tem- 
perature of the hot comptonizing electrons (kTe) and by an 
asymptotic power law index (F) related to kTe and to the 
medium optical depth r through the relation: 



9 



+ 



1 



0-(i + i) 



1/2 



(2) 



(see e.g. lLightman fc Zdziarskilll987l '). As no high energy cut- 
off appears in the XMM-Newton spectrum, we fix kTe to 
an arbitrary value of 100 keV. Such a choice does not af- 
fect significantly the results obtained on the XMM-Newton 
dataset alone. Interstellar absorption is treated with the 
phabs model, the a bunda nces are fixed at the values of 
lAnders fc Grevessd (Il989|). and cross s ection s are taken 
from iBalucinska-Church fc M cCammorJ (|l992l ) , and modi- 



fications after Yan et al 



([1998). Data show no significant 
deviation from solar abundances. The nthcomp component 
alone, poorly reproduces the observed spectrum (Xred = 
1.68 for 1604 d.o.f.). Residuals around 1.8 and 2.2 keV are 
evident. As there is no sign of similar features in the RGS 
spectra, we interpret them in terms of an incorrect calibra- 
tion of the instrumental Si and An edges that frequently 



ing that it is this component that provides the seed photons 
for Comptonization, we tie kTsoft to the temperature of the 
blackbody. The addition of these two thermal components 
definitely improves the fit, as the x^ decreases of 39 and 
486, for the addition of two and one degree of freedom, re- 
spectively. An absorption edge at 0.88 ± 0.01 keV is also 
clearly detected in the RGS data, with an absorption depth 
of 0.22 ± 0.06. We identify it as an absorption edge of O 
VIII {Eoviii = 0.871 keV). The edge is quite sharp, as an 
upper limit of 34 eV is found on its width, if a smeared edge 
(smedge) is used to model it. The best fitting parameters of 
this model (named as model A) are listed in the left column 
of Table [1 

The final reduced chi squared is anyway still large, 
Xred=l-26(1598 d.o.f.). Disentangling the temperature of the 
seed photons of the Comptonized component from the ob- 
served blackbody (model B, see Table [2)), significantly im- 
proves the fit (Ax^ = 232 for the addition of just one d.o.f.). 
Data thus favour a hotter and smaller region to produce the 
seed photons of the Comptonized component, than that of 
the observed blackbody. 

The presence of residuals around 6.6 keV (see Fig|6] 
where EPIC-pn residuals with respect to a simple power 
law are plotted in the bottom panel) suggests the presence 
of an iron Ka emission line. Adding a Gaussian centred at 



affect EPIC-pn spectra. An improvement of Ax = 151 is Epe — 6.65 ± 0.25 keV improves the x by 23 (for three 



obtained modelling these residuals with narrow Gaussian 
absorption features. 

As a soft excess is present, we add two thermal com- 
ponents to the model, accordingly to the results previously 
obtained mod elling XMM-Newton spectra of AMSP (GP05, 
IPapitto et al.i 2009). We model the softer component as disc 
emission (diskbb), and the hotter as a single temperature 
blackbody (bbodyrad). As the normalisation of the latter 
component is smaller than ~ 10 km, we interpret it as ther- 
mal emission arising from the neutron star surface. Assum- 



additional d.o.f.). The width of the feature {a > 0.56 keV) 
indicates it could be produced from reflection of the NS hard 
emission on the geometrically thin, optically thick, accretion 
disc. In this context, the feature is broadened by the rela- 
tivistic motion of the reflecting plasma in the inner parts of 
the accretion disc, where the space-time is bent by the gravi- 
tational influence of the compact object . We use a disklin e 
model to account for relativistic effects (JFabian et al.lll989l l , 
obtaining an improvement of Ax^ = 13 for 3 d.o.f. added. 
The diskline model describes the line shape in terms of the 
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size of the illuminated disc (that is, of its inner radius, Rtn, 
and of its outer radius, Rout), of the index of the radial de- 
pendence of the line emissivity, /3, and of the inclination of 
the system, i. As Rout is poorly constrained by the available 
statistics, we fix it to an arbitrary value of 10^ Kg (where 
Rg = GMi/(? is the NS gravitational radius, and Mi the 
mass of the compact object), of the order of the circular- 
ization radius for a system like IGR J17511-3057 (see e.g. 
iFrank et aLlbOOa ). The parameters we obtain are listed in 
the rightmost column of Table [2] (model C), and are compat- 
ible with the assumption that the line is emitted by plasma 
illuminated by the NS emission, and rotating in the accre- 
tion disc. The energy of the transition indicates mildly to 
highly ionised iron. As we show in the next section, these 
indications are supported also by the modelling of disc re- 
flection on a broader energy range. The feature stands at 
3.5cr above the continuum, and the F probability that the 
X^ improvement when switching from model B to model C 
is due to chance is < 10~^. However, the use of an F-Test to 
test the presenc e of additional compon ents has been strongly 
discouraged bv lProtassov et al.l (|2002l . but see also IStewarj 
|2009). We have then simulated 100 spectra starting from the 
best fit parameters of model B, and fitted them using both 
model B and C. In none of the cases we have noticed a x^ 
improvement, or an F statistics value, larger than the one 
we obtain from modelling real data. We therefore exclude at 
more than 99 per cent confidence level that the y^ improve- 
ment obtained with model C with respect to model B is due 
to chance. 

Model C is the best model we have found to fit the 
XMM-Newton dataset. However the x^ corresponding to 
this model is relatively large and the probability that we 
observe by chance a x^ equal or larger than the value we 
obtain, if the model is correct, is 0.5 per cent. This makes 
the model only barely acceptable. The average scatter of 
data points with respect to this model is anyway compatible 
with the accuracy of the effective area calibration of the 
EPIC-pn, while operated in fast modes (quoted to be 
better than 5% from an analysis of the 1.5-3 keV band, see 
fhttp://xmm2.esac.esa.int/docs/documents/C AL-TN-0083.pdf) 
Since the distribution of residuals does not show any sys- 
tematic trend, but are instead randomly distributed around 
what is predicted by our best fit model, we conclude 
that the x^ "^^ obtain is affected by uncertainties in the 
instrument calibration, and/or by the possible presence of 
unresolved and unfitted features. We therefore retain model 
C as a reliable description of IGR J17511-3057 spectrum in 
the considered bandwidth. 



3.3 A simultaneous RXTE-XMM-Newton spectrum 

The spectrum of IGR J17511-3057 is dominated by a power- 
law like component which we have interpreted in terms of 
Comptonization, with a cut-off temperature beyond the en- 
ergy band covered by the EPIC-pn. Moreover, the detection 
of a broadened Ka iron emission feature indicates how disc 
reflection may be important. To better assess these issues, 
we take advantage of the observations performed by RXTE 
(3-200 keV) overlapping with the XMM-Newton pointing. 
We thus added to the XMM-Newton data presented in the 
previous section, the PCA (3-50 keV) and the HEXTE (35- 
200 keV) spectra. A systematic error of 0.5% has been added 
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Figure 6. Residuals of the RXTE-PCA (top panel) and of the 
EPIC-pn (bottom panel ) with respect to an absorbed power law 
model. The RXTE-PCA points are subtracted from an additional 
background representing the emission of the Galactic ridge. 



to the PCA points according to the guidelines stated in the 
description of the latest PCA response matrix generator q 

In order to check the inter calibration between the 
EPIC-pn and RXTE-PCA spectra, we simultaneously fit 
these two spectra with an absorbed power law. The RXTE- 
PCA spectrum shows large swings up to ~ lOcr with respect 
to EPIC-pn points, with a clear soft excess below 5 keV and 
a much more prominent emission feature at ~ 6.6 keV. In 
order to investigate if such a discrepancy could be due to the 
contamination of the Galactic ridge emission, we subtracted 
from PCA data the Galactic ridge spectrum (see Sec l2.2[I H A 
great improvement (Ax^ ~ 376.33 for 710 d.o.f.) is achieved 
with respect to the unsubtracted data. Most importantly, 
the residuals of the EPIC-pn and the PCA points are dis- 
tributed in the same way after the subtraction of the Galac- 
tic ridge contribution to the PCA spectrum (see Fig. |6]). 
Such a result gives us confidence about the reliability of 
the combined XMM-Newton-RXTE spectral modelling we 
describe in the following. The relative normalisations be- 
tween the various instruments and the EPIC-pn are left free. 
The same normalisations quoted in the previous section have 
been found for the RGS, while PCA/PN and HEXTE/PN 
vary in the ranges 1.24-1.25 and 0.6-0.8, respectively, de- 
pending on the particular model used. We note that the val- 
ues taken by the normalisation factors for different spectral 
models are all compatible with each other within the respec- 
tive errors (that are typically 0.006 and 0.1 for PCA/PN and 
HEXTE/PN, respectively). 

As expected, the analysis of the combined datasets gives 
similar results than those presented in the previous section 
for what concerns the low energy part of the spectrum. Two 
soft components are requested to model the observed data, 
in addition to an a = 1.56 ± 0.04 power law, cut off at an 



^ http://www.universe.nasa.gov/xrays/programs/rxte/pca/doc/ 
/rmf/pcarmf- 11.7/ 

^ No additional background was considered for HEXTE data, as 
the Galactic ridge emission above 15 keV is minimal. 
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Table 3. Best-fitting parameters of the 0.5—200 keV spectrum composed of the RGS (0.5—2.0 keV), 

EPIC-pn (1.4-11 keV), PGA (3.0-50.0 keV) and HEXTE (35-200 keV) spectra. The continuum mod- 
els are defined as follows: (i) Model B, phabsxedge(diskbb+bbodyrad +nthcomp); (ii) Model B-Refi, 
phabsxedgeCdiskbb+bbodyrad+nthcomp+diskline+rdblurXpexriv); (iii) Model PS, phabsxedge(diskbb+bbodyrad+compps); 
(iv) Model PS-Refi, phabsxedge(diskbb+bbodyrad+compps+diskline+rdblurXpexriv). All the fiuxes calculated are unabsorbed. 



Model 


B 


B-Rofi 


PS 


PS-Rofi 


nH (10^2 cm'2) 


1 n4+°-'^-* 

-1-U4_0.05 


n Qs+O-04 


1 on+°''3 

J^-UU_Q Q4 


0-95t'ofs 


T O VIII 


n 22+0°^ 


0.14 ±0.08 


0.22 ±0.07 


0.19 ±0.07 


kT,„ (keV) 


27+°-''^ 

'-'•^'-0.02 


r, 00 + 0.05 

^■■'^-0.02 


r. O1+0.03 
'^■'^^-0.02 


0.36 ± 0.02 


Ri„v'cosi (dg km) 


33 ±9 


19 ±8 


24ll 


15^2 


FcLc(xlO"^erg/cm2/s:i^ 


0.20 


0.13 


0.17 


0.14 


Wbb (keV) 


0.53«:°? 


0.60 ±0.04 


0.60 ±0.02 


n fi4+ooi 


R-BB (^8 km) 


^■'-0.7 


6.8 ±0.6 


7.6 ±0.5 


6 3+°-3 
"■■^-0.1 


F°^|-««(xlO-9erg/cmVs) 


0.10 ±0.02 


0.10 ±0.02 


0.10 ±0.01 


0.11 ±0.01 


r 


1-^0-0.02 


1.78 ±0.02 


— 





kT„(keV) 


lOO(fixed) 


lOO(fixed) 


75+" 


51+^ 


rb 


0.6-1.9 


1.1-2.0 


1 10+0.08 
J^-t»_0.07 


1-34«:S? 


islab{°) 


— 


— 


65.0(/ixed) 


57.0+11 


kT^„j-t 


1 m+004 


1.1 ±0.1 


1.47 ±0.06 


-, oy + O.Ol 

^"^'0.02 


KsojMi km)^ 


4.2-5.3 


3.2-4.7 


q O + 0.3 

0.0_g 2 


o .+0.1 

■^■^-0.4 


Fg'-''-l''')(xl0-9erg/cm2/s) 


1.00 ±0.04 


1.04 ±0.07 


1.30 ±0.08 


1.00 ±0.05 


log 5 


— 


3.4 ±0.5 


— 


q n + 0.4 


n/2-K 


— 


0.18 ±0.05 


— 


14+01° 
'-'•^^-0.01 


F(,"/^7^^°'(xlO-Jerg/cmVs) 


— 


0.07 ±0.02 


— 


0.07 ±0.03 


EFe (keV) 


— 


6.8i«:? 


— 


6.82«:?? 


P 


— 


-Q+_l 


— 


-5.5l°1 


R„ {GM/c^) 


— 


22+_l 


— 


27tl 


i(°) 


— 


4ltl 


— 


48+?o 


EW (eV) 


— 


35 ±1 


— 


40.9 ±0.7 


Fjr^rCxlO-'-g/cmVs) 


1.19 ±0.04 


1.28 ±0.08 


1.52 ±0.08 


1.28 ±0.06 


1993.1/1721 


1838.1/1713 


1912.7/1720 


1847.7/1712 



'^ No error estimate is given as the unabsorbed disk fiux is evaluated down to 0.01 kcV, that is below the lower energy bound 
covered by data. 



for models B and B-Refi the range of optical depths has been evaluated from the best-fitting values of F obtained fixing kTe 
values in the range 50-200 keV, and using Eq. l[2}. 
■^ for models B and B-refi it has been evaluated assuming Fc = {e'^ — l)i^so/t (see EqlSj. 



at 



energy E — 51_g keV. Using nthcomp to model the hard 
emission we find that, if the observed blackbody provides 
the seed photons upscattered in the Comptonizing medium 
{kTsB = kTsoft = 0.34t°°^ keV), chi squared is large 
iXred = 1-26 for 1722 d.o.f.). Moreover, a region of size larger 
than the one measured (Rbb ~ 8.6 ±2.7 dg km, where dg is 
the distance to the source in units of 8 kpc) is needed to pro- 
vide enough photons for the observed hard component. Such 
a radius is evaluated considering that the flux escaping the 
Comptonizing medium is F' = j4 x Fsoft, where Fsoft is the 
seed flux, and the Compton amplification factor is ^ = e'', 
with y = {'ikTs/nieC?) x max (r,r2), the Compton parame- 
teijj. As the nthcomp model only evaluates the fiux added by 



^ The better this approximation holds the lowest the average 
energy of scatter ed photon is with respect t o the electron tem- 
perature, see e.g. iRvbicki fc LightmanI 1119791 ). 



the Compton process, Fc, the energy conservation between 
the two phases can be rephrased as Fc/Fsojt — e" — \. 
Expressing the seed fiux as Fsoft ~ crT^oftiRsoft/dY' , the 
radius of the area providing the seed photons has to be: 



Rsoft = 2.4 X 10"d, 



^Fc/jev ~ 1) 
{kTsoft^ 



km. 



(3) 



where kTsoft is expressed in keV. As using nthcomp no sig- 
nificant cut-off is found, we are forced to measure the best- 
fitting values of F fixing the electron temperature at values 
in the range 50-200 keV, and then evaluate r and y using 
Eq.Q. Estimating the flux in the nthcomp component in the 
0.5-150 keV energy range, the radius we obtain from Eq.© 
is Rsoft ^ 30 dg km, which is at least a factor ~ 3 larger 
than the normalisation of the observed blackbody compo- 
nent. We thus conclude that the observed blackbody is too 
cold to be the source of the seed photons. 
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Figure 7. Unfolded 0.5-200 keV spectrum of IGR J17511-3057 convolved with tiie model B+refl (left) and PS+Refl (right; see 
caption of Table |3] for a description), together with residuals in units of sigmas in the respective bottom panels. The overall 
model (thick solid line), the diskbb (dashed line), the blackbody (solid thin line), the Comptonized spectrum (dotted line), and 
the reflected spectrum (continuum and diskline, dash-dotted line) arc also overplottcd. Data have been re-binned for graphical 
purposes. 



Letting VTsoft 7^ kTss improves the modelling by 
Ax^ = 186 for the addition of one d.o.f. (see the leftmost 
column of Table [3] for the best fitting parameters of this 
model, named as B). As the electron temperature is not 
constrained by the model, we present the best fitting pa- 
rameters for kTe;= 100 keV. Using Eq.((3l) to evaluate the 
radius of the unseen thermal component providing the seed 
photons up-scattered in the Comptonizing medium, we find 
Rsoft —4-5 dg km. This radius is compatible with the ex- 
pected size of the hotspots on the NS surface. 

The presence of an iron Kq emission line in the XMM- 
Newton spectrum indicates a significant presence of reflec- 
tion. To model Compton reflec tion from a n ionised disc we 
use the pexriv model (Magdz iarz fc Zdz iarski 1995). The 
shape of the power law that describes the illuminating flux 
is fixed using the index F of the nthcomp component. The 
disc temperature is kept fixed at a value of 1 x 10® K, the 
disc inclination at the value indicated by the iron line mod- 
elling, and abundances to solar values. The relativity effects 
expected for a Keplerian accretion disc rotating around a 
compact objec t are taken into a ccount with the smearing 
kernel rdblur IjFabian et al.lll989l ) . To include also the most 
important bound-bound transition expected (Fe Kq) in the 
EPIC-pn range, we add a diskline with a line energy free 
to vary in the 6.4-6.97 keV interval. The outer radius of 
the illuminated disc is fixed at 10^ Kg, as the fit is rather 
insensitive for variations of this parameter in a wide range 
of values. The addition of a reflection component is highly 



significant, as the x^ improves by 155 for the addition of 
9 parameters (see the model BRefl in Table [3] for the best 
fitting parameters). 

As the Comptonization model nthcomp is not able to 
constrain the electron temperature, the assumption we have 
made (kTej^lOO keV) has potentially an impact on the 
amount of reflection needed to model the spectrum. The 
value of kTej also influences the determination of the seed- 
photon region (see Eq. [HJ . The results obtained considering 
both the 0.5-11 keV and the 0.5-200 keV energy bands in- 
dicate that the seed photons come from a region which is 
hotter and smaller than the blackbody, and not directly ob- 
served. It is then desirable to include in the model those seed 
photons that are not scattered in the r ~ 1 Comptonizing 
medium, as nthcomp only accounts for up-scattered photons. 
For all these reason, as well as to check the dependence of the 
measured parameters on the Compto nization model chosen, 
we substitute nthcomp with compps 1 Pout anen &: SvenssonI 
1996). This model has in fact a number of advantages with 
respect to the nthcomp model: (i) it enforces the energy bal- 
ance between the Comptonizing plasma and the region that 
provides the soft photons; (ii) it evaluates the Comptonized 
spectrum by solving numerically the radiative transfer equa- 
tion, for different scattering orders, and for different geome- 
tries of the comptonizing medium. It thus includes those 
seed photons that go through the hot medium unscattered. 
We consider in the following a slab geometry; (iii) it uses 
the pexriv kernel to model disc reflection of the hard pho- 
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tons only. The reflected spectrum thus imphcitly takes into 
account the decrease of photons incident on the disc at low 
( ^ 2 keV) energies, while using the model B-Refl the inci- 
dent spectrum is approximated as a power law extending to 
low energies. 

When reflection is not included in the model we obtain 
the results listed in Table |3] (Model PS). The continuum pa- 
rameters do not vary much when using this Comptonization 
model. In particular, the size of the region that provides the 
seed photons confirms the results previously obtained. Us- 
ing this model, we are also able to constrain the electron 
temperature of the hot medium, in agreement with the cut 
off energy found using a simple power law to model the hard 
emission. 

The property (ii) quoted above introduces a dependence 
of the emergent spectrum on the angle isiah between the slab 
normal and the line of sight. Keeping isiab as a free parame- 
ter in spectral fitting yields values ^ 80° , that would imply 
the presence of eclipses, which are instead not observed. As 
we see in the following, smaller values are obtained when re- 
fiection is included into the model. To calculate the spectrum 
presented in Table [3] we thus fix isiab ~ 65° . The tendency 
of the compps model to find large values of isiab when re- 
flection is not included, can be tentatively interpreted as an 
indirect indication of how reflection is significant to model 
the broad band spectrum. As a matter of fact, without re- 
fiection, a larger fiux has to be accounted for by the Comp- 
tonized component. As GP05 pointed out, at moderate op- 
tical depths like the ones we detect in IGR J17511-3057 
(r ~ 1), the Comptonization spectrum depends on isiab 
mainly in terms of the amount of unscattered photons ob- 
served, oc exp{—T I cos isiab)- FoT these optical depths the an- 
gular dependence of the scattered photons is indeed rather 
fiat (see the curves relativ e to angles ^ 70° in Fig. 4 of 
ISunvaev fc Titarchuklll985l ). As the spectral shape (hence 
kTsoft) is well constrained by data at low energies, it seems 
that compps gives account of the large fiux at high energies 
assuming that we see the slab at very large angles. In this 
way only a very small fraction of the unscattered photons 
would be observed and the total amount of seed photons 
would be as large as requested to ensure energy balance be- 
tween the hot and the cold phase. Besides the significant 
improvement of the model when refiection is included (see 
rightmost column of Table O model PS-Refi), the fact that 
such an overestimate of isiab is not needed when the addi- 
tion of refiection decreases the fiux at high energies owing 
to Comptonization, definitely supports its presence. 

The parameters of the iron line we find with model PS- 
Refi are all compatible with those found when Comptoniza- 
tion is described by nthcomp. Also the ionisation state and 
the amplitude of refiection indicated by the refiection model 
are consistent within the errors. The models B-refi and PS- 
refi are plotted in Fig[7]together with residuals. The addition 
of a refiection component and of an iron line to models B and 
PS gives values of the F-statistics > 7. In order to test if such 
an improvement may be due to statistical fiuctuations, we 
have repeated the procedure described in Sec l3.3l simulating 
100 fake spectra using the best fit parameters of the mod- 
els without reflection features. As we have never obtained 
an F-statistics value similar to that quoted above, when the 
refiection features are added to the model, we conclude at 
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Figure 8. Countrate recorded by the EPIC-pn in the 0.3-12 keV 
interval (a), apparent radius (b) and temperature (c) of the best 
fitting absorbed blackbody, fractional amplitude (d) and phase of 
the pulse first harmonic (e), during the first burst observed by the 
EPIC-pn. Time is given in seconds since the burst onset (MJD 
55095.05025). The values taken by the countrate, the radius and 
the temperature of the blackbody component of model PS-Refl, 
during the persistent emission are also over-plotted as dashed lines 
for comparison. 



more than 99% confidence level that the improvement ob- 
tained with their addition is not due to counting statistics. 
The best fit models presented for the combined RXTE- 
XMM-Newton dataset are only barely acceptable (null hy- 
pothesis probability of 1.8% and 1.2% for the models B-Refi 
and PS-Refi, respectively). Similarly to what has been noted 
in the previous section, the absence of systematic trends 
in the residuals leads us to conclude that the reduced chi 
squared is increased by uncertainties in the instruments cal- 
ibration and/or by the possible presence of unresolved and 
unfitted features. The fact that models obtained from a fit of 
five different instruments give indeed fairly satisfactory re- 
sults seems anyway worthwhile to note, and gives confidence 
about their reliability. 



4 THE BURST EMISSION 

IGR J17511-3057 exhibited two type I X-ray bursts during 
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the XMM-Newton pointing, sharing similar observational 
properties (see Table [4]). While the rise lasts for ~ 2 s, the 
decay follows an exponential decay {Xred =0.96 and 1.20, 
respectively, over 216 d.o.f.; see the top panel of Fig[8] for 
the light curve of the first burst). In order to study their 
temporal evolution, we extracted spectra over time inter- 
vals or variable length (from Is at the beginning to 8s at 
the end), subtracting as a background the persistent spec- 
trum analysed in the previous section. Pile-up is not of con- 
cern, as it affects energies ^ 8 keV for at most few per 
cent. We have anyway checked that the spectrum extracted 
without the brightest CCD column gives spectral model pa- 
rameters compatible with those quoted hereafter. We suc- 
cessfully model the burst emission with an absorbed black 
body. The absorption column has been varied in the range 
indicated by the various models we have used to model the 
persistent spectrum. The evolution of the temperature and 
of the emission radius of the first burst, are plotted in pan- 
els (b) and (c) of Fig[8l Results from the second burst are 
qualitatively similar. The burst temperature follows the ex- 
ponential decay of the burst flux; the apparent emission ra- 
dius initially increases but soon reaches an asymptotic value, 
suggesting that the critical threshold for photospheric ra- 
dius expansion was not reached. In the time interval 6-30 s, 
the apparent radius reaches approximately a constant value, 
Rapp = (7.0±0.3)dg km. We restrict to the 6-30 s interval to 
estimate this radius as subsequently, the burst flux becomes 
comparable with the persistent flux, and systematic errors 
could arise due to the subtra ction of the persistent emission 
jvan Paradiis fc Lewin|[l986l '). 

In panels (d) and (e) of Fig[5] the amplitude and the 
phase of the pulse profile computed on the first harmonic, 
are also plotted. Data have been preliminary corrected with 
the best orbital solution and folded in 8 phase bins around 
the spin frequency derived in Sec l3.1l During the first 4s 
of the burst, the amplitude decreases to a value compati- 
ble with zero. Though, as the telemetry limit of the EPIC- 
pn is trespassed at those countrates, the disappearance of 
pulsations may also be due to a dead time effect. We then 
discard pulse data referring to the first few seconds after 
the burst onset. However, pulsations are soon recovered and 
observed throughout the decaying part of the burst. The am- 
plitude is roughly comparable with that of persistent emis- 
sion, and also the phase is stable within ~0. 1-0.2 cycle. This 
behaviour suggests that a mechanism similar to the one orig- 
inating the persistent pulsations is at work during the type-I 
X-ray burst (see Sec l5.4p . 



5 DISCUSSION 

5.1 The X-ray pulsations 

X-ray coherent pulsations at a frequency v — 244.8339512(1) 
Hz are observed throughout the observation performed by 
XMM-Newton. No spin frequency evolution is detected at 
a 3 (T upper limit of 1 x 10~^^ Hz s~^. Such limit sim- 
ply refiects the shortness of the time interval covered by 
the XMM-Newton observation. As a matter of fact, if a NS 
spinning at the frequency of IGR J17511-3057 simply gains 
the specific angular momentum of the accreted matter, the 
expected spin frequency derivative cannot be larger than 



Table 4. Parameters of the type I X-ray bursts observed by 
XMM-Newton. Errors in parentheses are evaluated at a 90% con- 
fidence level. 



Burst I 



Burst II 



T,tart (MJD) 


55095.05025 


55095.52363 


r (s) 


11.4(2) 


11.5(2) 


poo (a) 

peak ^ 


4.2 ±1.1 


4.3 ±1.1 


_poo (b) 


4.8 ±1.1 


4.9 ±1.2 



^, evaluated from 
the best fitting blackbody spectrum at the peak of the outburst. 
Burst fluence in units of 10"'' erg cm~^ 



2.4 X lO-^^IZ^LsrRiom^^/! Hz/s. 



(4) 



Here I45 is the NS moment of inertia in units of 10 g cm , 
Til, 1.4 is the NS mass in units of 1.4 Mq, Riq is the NS 
radius in units of 10 km, and L37 is the bolometric X-ray 
luminosity in units of 10''^ erg s~^, the dependence on which 
have been introduced using the relation Lx ~ GMM/Rns ■ 
Thus, considering the estimate of Lx we give from spectral 
modelling (~ 1.6 x 10'^^ d| erg s~^) the maximum expected 
spin frequency derivative is nearly two order of magnitude 
lower than the loose upper limit we could set, given the 
length of the XMM-Newton observation. 

The pulsed fraction is among the largest ever observed 
from an AMSP (14.4(3) % RMS). Typical observed values 
are between 2 and 8 %, with sporadic increases up to 
^ 10% especially at the end o f the outburst episodes 
( Wiinands fc van der Kli^ Il998l: iMarkwardt et al.l I2OO2I: 



2005T^G allowav et al 



1 T 



i 



Altamirano et al. 



20051. 



2007: 



Galloway et al.1 |2002|: ICampana et al. l2003l: IWatts et all 



Casella et al.l I2OO8I : 



Patruno et al.1 bOld ). 



while 
amplitude up 



Patruno et al.l ([20091) detected a pulse 
to 19% from XTE J1807-204. If the misahgnment angle 
between the magnetic dipole and the spin axis is small 
{9 ^ 20°), the observed value of t he pulsed fraction suggests 
that the inclination is ^ 45° (IPoutanen fc Beloborodovl 
120061 ). Nevertheless, the degree of anisotropy of the emitted 
light and the assumed spot shape may significantly influ- 
ence the pulse amplitude. A detailed analysis of the pulse 
profile is mandatory to derive firm constraints on the spot 
geometry, and as it is beyond the scope of this paper, we 
defer it to another work. It is worth to note, however, that 
the presence of a hump in antiphase with respect to the 
global maximum possibly suggests that we see, at least for 
a fraction of rota tional phases, the antipodal spot (see e.g. 
iLeahv et aLll2009l ). 

The pulse profile of IGR J17511-3057 shows strong 
spectral variability, similaxly to what has been observed from 
the AMS Ps, XTE J1751-305 and SAX J1808.4-3658, by 
GP05 and lPatruno et al.l ((20091), respectively. The decrease 
of the pulsed fraction at low energies can be completely 
ascribed to to the growing influence of the accretion disc 
emission. To see this, we plot in FigOJthe fractional contri- 
bution of the bbodyrad and the compps component (Model 
PS-refl), to the overall EPIC-pn countrate, f(BB+PS) (top 
curve, right scale). The flrst harmonic fractional amplitude 
is also plotted (bottom curve, right scale, same as Fig|4]). 
The similarity between the shape of the two curves is strik- 
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Figure 9. Fractional contribution of the bbodyrad and compps 
components to the EPIC-pn countratc (top curve, left scale), and 
first harmonic fractional amplitude (bottom curve, right scale). 



ing. The correlation can also be quantitatively expressed, as 
the ratio between the amplitude in the range 1-2 keV and 
the maximum amplitude (ylf "^ = 0.232(7), E=9-ll keV) is 
0.71(2), almost equal to /(bb+ps) in tlie same energy band 
(0.675(3)). The energy dependence of the pulsed fraction 
therefore represents a compelling evidence that the bulk of 
the emitted spectrum is somewhat related to the magnetic 
caps, where the pulsations are formed. The similarities with 
the pulse amplitude energy dependence of other AMSPs (see 
above), further strengthens such a conclusion. 

Soft phase lags of ~ 230/is have also been observed. Soft 
lags of a similar amount (~ 20 0^ts) were observ ed from the 
AMSP SAX J1808.4-3658 bv lCuiet all (|l998), who inter- 
preted them in terms of Compton downscattering of intrin- 
sically hard photons in a r ~ 10 cloud, which is however not 
observed in the X-ray spectrum. It was also proposed that 
soft lags could be due to the Doppler energy shifts intro- 
duced by the NS fast rotation (lFordll200Gl : [Weinberg et al.1 
[200l|), or due to the different angular distribution of the flux 
emitted by the blackbody and the Comptonized emission 
l|Gierliiiski et al.ll2002l : [Poutanen fc Gierliiiskill2003l '). As the 
lags we observe from IGR J17511-3057 are of similar am- 
plitude than those of SAX J1808.4-3658 (which spins 1.7 
times faster than this source), Doppler effects are unlikely 
the only physical reason behind phase lags. 



5.2 The companion star 

IGR J17511-3057 is one of the AMSP with the longest or- 
bital period ever observed. From the measured mass func- 
tion, f{m2;mi,i) = (m2 sini)^/(mi -|- m2)^ (see Table [T|, 
it is possible to derive constraints on the mass of the com- 
panion star, 1712. Here f,mi and 1712 are expressed in solar 
masses. The absolute minimum on 7712 is reached consid- 
ering an inclination of 90°, which translates in m™"" ~ 
0.05 -I- 0.06mi for 1 ^ mi ^ 2. The absence of eclipses lim- 
its the range of viable inclinations to i < 90° — 7, where 
7 — arctan (i?2/a) is the angle subtended by the companion 
star as seen from the NS, R2 is the radius of the compan- 




0.05 



0.3 0.35 

"2(Msun) 

Figure 10. Size of the companion Roche Lobe (top curv e) and 
ZAMS mass radius relation as from Chabrier & Baraff^ ||2000|, 
bottom curve). The values taken by an assumed Roche Lobe fill- 
ing companion star, for a set of inclinations and mi = 1.4, are also 
indicated. The left shaded region is excluded from the absence of 
eclipses, while the right shaded one from the assumption that the 
companion is not larger than its Roche Lobe. The minimum mass 
for i < 70° is also plotted with a vertical dashed line. 



ion star and a is the binary separation. In order for mass 
transfer to proceed the radius of the companion must be 
comparable with the radius of the Roche Lobe, R2 ^ Rl2, 
which is estimated through the relation given bv lPaczvhskil 
(|l97ll ): 



Rl2 — 0.462a[m2/(mi -f- m2 



,1/3 



(5) 



To evaluate the maximum inclination compatible 
with the absence of eclipses, we calculate 7 for R2 = 
RL2{m2 ), obtaining i ^ 78°. The absence of dips in the 
X-ray Ught curve further decreases the range of possible in- 
clinations, so that we consider i ^ 70° only, slightly increas- 
ing the minimum companion mass to, e.g., m2 ^ 0.15, for 
mi = 1.4. A main sequence companion star is clearly in- 
dicated. To see this we use the third Kepler law in Eq.© 
to relate the Roche Lobe radius to the mean density of the 
companion star at the observed orbital period: 



rL2{m2) = Rl2/Rq = 0.536 m. 



1/3 



(6) 



This relation is plotted in Fig llOl together with the 
Zero Age Main Sequence (ZAMS) mass-radius relation of 
IChabrier fc Barafld (|2000l '). rf-^(m2). The two curves cross 

— 0.41. Heavier companion 



for m 



(max) 



0.44 and r. 



(max) 



stars can be safely excluded as the Roche Lobe would be 
overfilled, and mass transfer unstable. Assuming mi — 1.4, 
this solution corresponds to i^'""') = 20.5°. 

The reflection continuum and the iron Ka line we have 
modelled with a diskline indicate an inclination in the 
range 38-68° (see Table |3] and Sec l5.3p . This interval corre- 
sponds to 7712 =0.15-0.23 for mi = 1.4. Slightly larger values 
are found for heavier NS. The inclination estimate we have 
drawn from spectral fitting thus indicates that the compan- 
ion is slightly bloated with respect to its ZAMS thermal 
equilibrium radius. A possible mechanism to drive a com- 
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panion star out of thermal equilibrium is irradiation by the 
compact object, but t o what exte nt it is important is still 
matter of debate (see iRitteii 120081 . and references therein). 
However, irradiation is not strictly needed to account for 
the values we observe from IGR J17511-3057, as CV-like 
evolution can explain it if the companion star is slightly 
evolved at the onset of mass transfer ( but still before the 
MS t urn off), or heavier than 1 M©. As lPvlvser fc Savoniid 
l|l989) shown, angular momentum losses driven by magnetic 
braking and gravitational radiation can still lead to a con- 
verging system (that is, a system that evolves decreasing its 
orbital period), even for slightly evolved, or heavy, compan- 
ions. In these cases, the companion radius will be systemat- 
ically larger than the ZAMS radius. To see this, we consider 
the ev olutionary tracks calculated by IPodsiadlowski et al.l 
l|2002f) . While an unevolved companion with initial mass of 
mj =1 crosses the observed orbital period of IGR J17511- 
3057 when m2 — 0.44, as its radius is frozen to the ZAMS 
value, this happens for m2 — 0.16 if it is slightly evolved, 
and thus larger, at the beginning of mass transfer. Similar 
results are obtained even for initially unevolved, but heav- 
ier, companion star (the evolutionary track for a m2 =1.4 
companion star passes through the period of IGR J17511- 
3057 when m2 = 0.15). Values of m2 in agreement with the 
inclination indicated by spectral modelling can be thus ob- 
tained, for peculiar conditions of the companion star at the 
onset of mass transfer. 



5.3 The X-ray spectrum 



The correlation between the energy dependence of the 
pulsed fraction and the spectral decomposition we have em- 
ployed strongly suggests that the Comptonizing medium 



surrounds the hot spots o n the NS surface. Giorliiiski et al 



(2002), Pouta nen fc Gier lihski (2003), GP05. iFalanga et al 



We presented in Sec l3.3l and 13.31 a detailed spectral analy- 
sis of the 0.5-200 keV X-ray emission of the AMSP IGR 
J17511-3057, using two different Comptonization models, 
and checking the presence reflection. The spectrum of IGR 
J17511-3057 is modelled by four components, (i) the accre- 
tion disc emission peaking at ~ 0.3 keV; (ii) a ~ 0.6 keV 
blackbody of apparent radius ~ 7 dg km; (iii) Comptoniza- 
tion from a hot {kT^ — 50 keV) medium of moderate optical 
depth (r ~ 1.3), of the thermal (~ 1.3 keV) photons pro- 
vided by a ~ 3.5 dg km region (iv) Compton reflection and 
an iron Ka emission line, interpreted as coming from the ac- 
cretion disc illuminated by the hard radiation emitted from 
the NS surface. The shape of the best fitting model does not 
change much using different Comptonization models. As the 
model PS-refl better addresses a number of physical aspects 
of this source emission (see Sec l3.3|l we discuss the param- 
eters thus obtained. Compton disc reflection is significantly 
detected regardless of the particular model used to describe 
Comptonization, as well as an iron line at an energy compat- 
ible with Fe XXV-XXVI. Even if reflection has not always 
been observed, such a spectral decomposition have already 
proved s uccessful in the descr iption of the X-ray spectrum 
AMSP iGierlinski et all I2OO2I, GP05, iFalanga et all I2OO5I : 
iPapitto et al.ll2009l : [Patruno et al.l 120091 '). In particular, the 
presence of disc thermal emission has always been detected 
when XMM-Newton high resolution spectra extending down 
to ~1 keV were available. The similarities in the X-ray emis- 
sion of AMSP can be viewed as an indication of how the 
physical processes that produce the bulk of their spectra 
are similar, and probably related to the presence of a mag- 
netosphere in these systems. 



(|2005i) and iFalanga et al.l (|2005|) interpreted the similar hard 
components shown by other AMSPs as coming from Comp- 
tonization in a plasma heated in the accretion columns. The 
magnetic field coUimates the in-falling matter in columns 
of radii that are few tenths of the NS size. Even for accre- 
tion rates ~ 0.1 Msdd like those observed for these sources 
during outbursts, the local accretion rate rh can the refore at- 
tain the Eddington level IjBasko fc Sunvaevlll976r ). In such 
a case, the radiation pressure becomes comparable to the 
ram pressure of the plasma falling at supersonic velocities, 
and a shock may thus form, heating the electrons up to the 
observed large temperatures. 

The observed ~ 0.6 keV blackbody component is too 
cold to supply enough seed photons for the observed Comp- 
tonized spectrum. Disentangling the temperature of seed 
photons from that of the observed blackbody significantly 
improves the fit and ensure the energy balance between the 
hot plasma and the region that provides the seed photons. A 
similar result has already been obtained by GP05, modelling 
the X-ray spectrum of the AMSP, XTE J1751-305. The size 
of the region that provides the soft photons, K^aft ~ 3.4^g'4 
dg km, is compatible with the expected radius of an hot 
spot. The radius of the observed blackbody, R^b ~ 6.3J;q 'J 
dg km, is instead of the same order of the asymptotic value 
taken by the radius of the thermal burst emission. This com- 
ponent is therefore emitted by a larger fraction of the NS 
surface. 

The cooler thermal component can be safely attributed 
to the emission of an optically thick accretion disc. The con- 
tribution of this component at low energy is responsible for 
the decrease of the pulsed fraction. Moreover, the parame- 
ters we obtain with a diskbb model nicely fit the expecta- 
tions for an accretion disc around an AMSP. Considering the 
range of inclinations indicated by the refiection component 
(38-68°), the model diskbb evaluates an apparent inner disc 
radius Rin in the range 18-28 dg km. Such a radius is esti- 
mated approximating the disc temperature profile with its 
asymptotic behaviour T(r) ex: r~^'* , without accounting for 
the decrease of the viscous torque at the inner boundary 
of the disc, and neglecting spectral hardening. To evaluate 
the imp ortance of these effec ts we consider the disc model 
diskpn ICierliiiski et al.ll 19991 ). The spectral shape found by 
diskbb is recovered after rescaling the inner disc radius Rin 
by a factor « 2.3, for mi = 1.4, d = 8kpc, i = 50° and 
/ = 1.7, where f is the ratio between the col our and ef- 
fective temperature (|Shimura fc Takaharall 19951 ). The inner 
radius indicated by disc emission modelling thus meets the 
request that an accretion disc around an X-ray pulsar is not 
truncated at radii much larger than the corotation radius. 
Re = [GM/A-K'^u'^f^^ = 42.8mi.4 km for IGR J17511-3057, 
where mi. 4 is the NS mass in units of 1.4 Mq . Also the ratio 
between the flux in the disc component and the total flux 
agrees with a truncation radius a few tens of km away from 
the NS star. From the virial theorem it is in fact expected 
that Fdisc/Ftot ^ i?*/2ri„. As we measure Fdisc/Ftot ^ 0.1, 
the fiux ratio predicts r^ ~ 5R„ , which substantially agrees 
with the measured value. 
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The presence of reflection features, such as an iron Ka 
emission line and Compton reflection, is significantly de- 
tected in spectral fitting. The iron line stands at '^ 4cr above 
the continuum and has an equivalent width of ~ 40 eV. 
If modelled with a Gaussian, its width can only be con- 
strained to be ^ 0.5 keV. As the spectrum of this source is 
dominated by Comptonization in a hot and optically thin 
medium, it is hard to find alternatives to disc reflection in 
order to explain such a width. In the inner regions of the ac- 
cretion disc the line is broadened and red-shifted by the ef- 
fects of Keplerian motion and by the gravitational influence 
of the nearby NS. Modelling the feature with a diskline wc 
obtain an estimate of the inner disc radius as ri„ = 27lg Rg, 
which translates in ri„ = 561^9 km for a 1.4 M© NS. Even if 
it is rather loose due to the limited counting statistics, also 
this estimate overlaps with the corotation radius, as it is ex- 
pected from accretion theories. The addition of the RXTE 
dataset to the XMM-Newton spectrum also allows the de- 
tection of the hump at --^ 30 keV expected from Compton 
reflection. The ionised state of the reflector is rather high, 
log^ = S.OIq'j. It is worth to note that the highly ionised 
surface in the refiector is in agreement with the transition 
energy of the broad iron line at Efe ~ 6.82lg]^j^ keV, that 
we identify as a Lya resonant transition of Fe XXV (rest- 
frame energy at 6.70 keV), that is li kely to be produced i n 
a photo-ionised plasma at log^ ~ 3 (|Kallman et al.ll2004h . 
Considering the values indicated by the reflection contin- 
uum and by the iron line modelling, we give a conservative 
estimate of the inclination in the range 38-68°. Also the 
OVIII edge we find at a high absorption depth (r =0.1-0.3) 
can be tentatively interpreted in terms of disc reflection. As 
its width is < 34 eV, the outer rings of the disc should be 
involved, in order to make the rotational broadening negli- 
gible. 



5.4 The X-ray bursts 

In Sec. |4] we have presented a spectral and temporal anal- 
ysis of the bursts obse rved by XMM -Newton. Considering 
the Swift observations, IBozzo et al.l 12010.) concluded that 
the observed recurrence times could be reconciled with the 
persistent flux if it is assumed that the bursts are ignited in 
a pure helium environment. 

This hypothesis can be checked estimating the local ac- 
cretion rate rn needed to produce the fluence of the sec- 
ond burst, which lags the first by At =40.9 ks. One has in 
fact that rn — y/At, where y is the column depth at which 
the burst is ignited. The val ue of y can be estim ated as 
y = Aird^T/QnucRt (see e.g. iGallowav et aLll2008l ). where 
!F = {1 + z)TF^^f. is the burst fiuence, and Q„uc is the en- 
ergy per nucleon released during the thermonuclear burning 
{Qnuc ~ 1.6 MeV fo r complete burning of He into iron group 
elements I Wallace fc Wooslcy 1981). The measured values 
(see TablelH) yield y = (2.6 ±0.6) x 10** dj R^q g cm'^, and 
m = (6.3 ± 1.4) X 10^ dj R'q g cm"^ s"^ ~ 0.08 R^QniEdd, 
where z — 1.31 is considered. This value of rh agrees with 
the bolometric luminosity of the persistent emission. Assum- 
ing that the blackbody and the Comptonized component 
arise from the vicinity of the NS surface, and therefore the 
observed fluxes have to be corrected for general relativity 
effects, we in fact estimate from model PS-refl the persistent 
luminosity as, L^""' = 4nd'^[{l + zf{FBB + Fps) + {FD^sc + 



FRefi)] = (1.6 ± 0.1) X 10^^ di erg s'^ ~ 0.09 mj-^ Ledd. 
The agreement between the persistent luminosity and the 
mass accretion rate needed to explain the bursts fluence and 
recurrence time, when it is assumed that all the hydrogen 
is steadily burnt between bursts, suggests that bursts are 
triggered in nearly pure helium environment. It is perhaps 
useful to clarify that the difference between the local accre- 
tion rate estimated by the burst properties, m ~ O.OSniEdd, 
and the rate needed for a shock to form in the accretion 
column, m ~ ruEdd, arises as the first is evaluated aver- 
aging the accretion rate over the whole NS surface, while 
the latter considers only the column size. Such a treatment 
seems appropriate as a magnetic field ^ 10^ G like the one 
expected in an AMSP should be able to confine accreted 
matter i n the polar caps up to a column density of ~ lO'^ 
g cm~'^ (JBrown fc Bildstenlll998r ). which is well below the 
value we measure. The accreted layer is therefore able to 
diffuse across the whole NS surface before the burst onset. 
As no photospheric radius expansion is observed, an up- 
per limit on the distance can be set as cf ^ LpRE/iTrF^^f., 
where the luminosity at which photospheric radius expan- 
sion is expected, Lpre = (3.79 =b 0.15) x 10** ** erg s~^, was 



empirically derived by iKuulkers et al.l (|2003r ) . Considering 
the peak fiux of the second burst and the associated uncer- 
tainty, we set an upper limit on the distance to the source of 
10 kpc at a 90% confidence level. The radius of the thermal 
emission attained during the burst decay can be identified 
with the NS radius, and used to put a reasonable lower 
limit on the distance to the source. The measured value 
Rapp = (7.0 ± 0.3)^8 km can be related to the effective ra- 
dius through the relation, R = Rappfci^ + z)~^ (see e.g. 
iLewin et al.lll993r ). The factor fc takes into account spec- 
tral hardeni n g, an d has been estimated as fc — 1.35 by 
iMadei et al.l (200J), for the case of a source that does not 
reach the Eddington level. Solving the previous relation for 
M > 1.4 Mq, and considering that a soft equation of state 
like the A predicts 7?* > 8.5 km and M < 1. 92 Mp for a NS 
spinn i ng at the rate of IGR J17511-3057 IjPandharipandel 
Il97ll : ICook et al.]|l994h . we obtain d> 6.5 kpc, unless a 
strange star is considered. For a lighter NS (M=1.57 M©) 
this limit increases to 6.9kpc. Given the estimated range of 
distances, it is possible to safely state that IGR J17511-3057 
belongs to the Galactic bulge. Also the compatibility of the 
measured equivalent hydrogen absorption column with that 
observed from XTE J1751-305 (whose distance was esti- 
mate d to be ^ 7 kpc. lMarkwardt et al.ll2002l : |Papitto et al.l 
[2003) supports such a conclusion. 

Pulsations are detected throughout the burst decay at 
a similar amplitude than that seen during the persistent 
emission, while during the first few seconds they disappear, 
probably because of the telemetry limitations of the EPIC- 
pn. Their phase is also compatible within 0.1-0.2 cycles 
with that of pre-burst pulsations, suggesting that magnetic 
channelling of accreted matter and the burst onset hap- 
pen at a similar location on the NS surface. This indication 
is strengthened by the analysis performed by Riggio et al. 
(2010, submitted) over RXTE data. They find in fact that, 
with the exception of the first two seconds since the burst 
onset, the burst oscillations are phase locked to pre-burst 
pulsations within 0.1 cycle. Phase locking between burst 
and non-bu rst oscillations of a n AMSP has already been 
observed bv lWatts et al.l (|2008r ). at a much larger accuracy 
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(~0.01 cycles) than that made available by the data pre- 
sented here. As they discuss (see also references therein), 
the most appealing interpretation is in terms of the effect 
a temperature gradient between the fuel impact point on 
the surface, and the rest of the surface, may have on the 
ignition conditions, or on its development. A temperature 
gradient during the persistent emission is indeed indicated 
by our spectral modelling, which requires the region feed- 
ing the accretion columns of soft photons to be hotter than 
the rest of the surface. A firm explanation of the correlation 
between pulsations observed during bursts and during the 
persistent emission is anyway still missing. 
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